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O ' Abstract. The steady-state of a generahzed coagulation-decoagulation model on a 

onc-dimcnsional lattice with reflecting boundaries is studied using a matrix-product 
approach. It is shown that the quadratic algebra of the model has a four-dimensional 
representation provided that some constraints on the microscopic reaction rates are 
fulfilled. The dynamics of a product shock measure with two shock fronts, generated 
by the Hamiltonian of this model, is also studied. It turns out that the shock fronts 
C^ ' move on the lattice as two simple random walkers which repel each other provided that 

ff , the same constraints on the microscopic reaction rates are satisfied. 
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1. Introduction 

The study of microscopic structure and dynamics of traveling wave solutions in multi- 
species one- dimensional stochastic systems has attracted people's attention in this field 
considerably in recent years [l]-[15]. For instance, the microscopic dynamics of shocks 
are studied for three families of single-species one-dimensional reaction-diffusion systems 
with open boundaries and nearest-neighbors interactions which include the Partially 
Asymmetric Simple Exclusion Process (PASEP), the Branching- Coalescing Random 
Walk (BCRW) and the Asymmetric Kawasaki- Glauber Process (AKGP) P]. It has been 
shown that in all three systems the time evolution of a product shock measure with a 
single shock front is equivalent to that of a simple random walker on a finite lattice with 
homogeneous hopping rates in the bulk and special reflection rates at the boundaries, 
provided that some constraints on the microscopic reaction rates are fulfilled. The 
steady-states of these three systems can be essentially written as a linear superposition 
of such product shock measures. On the other hand, the steady-states of these systems 
can be obtained using the matrix-product formulation [16] in which the steady-state 
weights are written in terms of the product of non-commuting operators which satisfy 
a quadratic algebra (for a recent review of this approach see pJJ). Surprisingly, it 
has been found that the conditions under which these operators have two-dimensional 
matrix representations are exactly those for a product shock measure to have a simple 
random walk dynamics in these systems [16] . 

It has been shown in [5J that the existence of a two-dimensional representation for 
the quadratic algebra of a multi-species one-dimensional reaction-diffusion system with 
open boundaries and nearest-neighbors interactions implies that the steady-state of the 
system can be written in terms of a linear superposition of product shock measures. 
It seems that the same is true for the systems whose quadratic algebras have higher- 
dimensional matrix representations. For instance, for the Totally Asymmetric Simple 
Exclusion Process with open boundaries it is known that the quadratic algebra has 
an infinite-dimensional matrix representation. It has been shown that this is associated 
with the fact that the steady-state of this system can be written as a linear superposition 
of product shock measures with infinite number of shocks [6]. 

The microscopic dynamics of shock fronts has also been studied in driven-diffusive 
systems with more than a single species of particles and also in these systems with 
next-nearest-neighbor interactions [T0]-[T5j. 

Despite of these efforts it is still generally unclear whether the existence of a finite- 
dimensional matrix representation of the quadratic algebra is related to the fact that 
the steady-state can be expressed as a superposition of product shock measures. 
In this paper we consider a generalized one-dimensional single-species coagulation- 
decoagulation model with reflecting boundaries as a new example. We believe that 
the study of this exactly solvable model provides us with another piece of evidence for 
the existence of such relation and definitely sheds more light on the unknown aspects 
of this problem. We will show that a product shock measure with two shock fronts 
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which have simple random walk dynamics can evolve in this system provided that the 
microscopic reaction rates satisfy some constraints. This will enable us to construct the 
steady-state of the system simply by considering a linear superposition of such measures; 
however, we will not follow this approach here. Instead we will show that under the same 
constraints the quadratic algebra of this model has a four- dimensional representation. 
Our paper is organized as follows: In section (2) we will explain the mathematical 
preliminaries and introduce the model. In section (3) we will study the dynamics 
of a product shock measure with two shock fronts and investigate the conditions (by 
imposing some constraints on the microscopic reaction rates of the model) under which 
it has a simple time evolution similar to that of two random walkers moving on a one- 
dimensional lattice while reflecting from the boundaries. The fourth section will be 
devoted to the investigation of whether a finite-dimensional representation exists for 
the quadratic algebra of the model under the same conditions. In the last section we 
will discuss the summery of results. 

2. The model 

We define \P{t)) as the probability vector of a Markovian interacting particle system at 
the time t. The time evolution of this vector is governed by the master equation which 
can be written as a Schrodinger equation in imaginary time as follows: 

j^\P{t)) = H\P{t)) (1) 

in which H is an stochastic Hamiltonian. The matrix elements of the Hamiltonian are 
the transition rates between different configurations. For a single-species system with 
nearest-neighbors interactions defined on a one-dimensional lattice of length L with 
reflecting boundaries the Hamiltonian H is of the following form: 

L-l 

H = z2 hk,k+i (2) 

fc=i 

in which: 

where X is a 2 x 2 identity matrix and /i is a 4 x 4 matrix for the bulk interactions. In 
the basis (00, 01, 10, 11) the Hamiltonian h can be written as: 

/ LOn UJi2 CJi3 UJi4 \ 



h 



UJ21 UJ22 "^23 ^24 

UJ31 UJ32 ^33 ^34 

\ UJ41 UJ42 ^43 ^44 / 

in which we have defined the basis vectors: 



(3) 



associated with the presence of a particle and a hole in each lattice site respectively. 
Requiring the conservation of probability, the sum of the elements in each column should 



Finite- dimensional representation of the quadratic algebra 4 

be zero; therefore, one has Ua = — Y^j^i ^ji for hj = 1; ' ' ' ; 4. 

In our generahzed coagulation-decoagulation model the non-zero reaction rates belong 

to the following reactions: 

+ A -* A + with rate cjga 

v4 + ^ + A with rate UJ23 

A + A ^ A + ^ with rate uju f^s 

A + A^iD + A with rate UJ24, 

(/) + A -^ A + A with rate a;42 

A + iJ) ^ A + A with rate uJi3 

in which A and stand for the presence of a particle and a hole in each lattice site 
respectively. This model has already been studied in related literatures in a couple of 
special cases. For instance consider the case: 

c^24 = ^23 = q'^ , ^34 = uj32 = q , c^42 = Ag , u;43 = Aq'^ . (6) 

The complete spectrum of the Hamiltonian in this case has been obtained exactly |18j . 
The steady-state of the system has also been studied using both the method of empty 
intervals pj] and the matrix-product approach [19] . It turns out that the phase diagram 
of the model in this case has two different phases: a low-density phase and a high-density 
phase which are separated by a coexistence line. The particle concentration on the lattice 
in both phases have exponential behaviors with three different length scales while on 
the coexistence line it changes linearly in the bulk of the lattice and grows exponentially 
near one of the boundaries. By using the matrix product approach it has been shown 
that for the special tuning of the parameters given in ([6]) the quadratic algebra of the 
model has a finite-dimensional matrix representation [19]. Since the system has three 
different length scales and that these length scales are determined by the eigenvalues of 
these matrices [18j, the authors in p^ have found that the minimum dimensionality of 
these matrices should be four. 

The most general case ([5]) has already been studied on a finite lattice with open 
boundaries in which the particles can enter into the system or leave it from the 
boundaries [4J and also on a lattice with periodic boundary conditions [8]. Note that 
the microscopic reaction rates in our generalized model are exactly those of the BCRW 
studied in P]. As we mentioned before, it is known that traveling wave solutions exist 
for the BCRW with open boundaries under some constraints on the microscopic reaction 
rates. In [16] the author has shown that in the open boundaries case the steady-state can 
be written as a matrix-product state with two-dimensional matrix representation given 
that the same constraints are satisfied. The authors in [9] have studied the dynamical 
phase transition in this system when considered on an infinite lattice. 
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Figure 1. A simple sketcli of a double-sliock measure. Tlie sliock positions are defined 
at the sites m and n. 



3. Dynamics of a product shock measure with two shock fronts 

In what follows we consider the most general reaction rates ([5]) and study the time 
evolution of a product shock measure with two shock fronts defined as: 

for < m < n — 1 and 1 < n < L + 1 

on a lattice of length L with refiecting boundaries. We define two auxiliary sites 
and L + l for convenience. In Figure 1 a simple sketch of a double-shock structure is 
given. We investigate the conditions under which the time evolution equations for ([71) 
are simply those of two random walkers at sites m and n moving on a finite lattice. It 
turns out that by imposing two constraints on the microscopic reaction rates one finds 
the appropriate answer. These constraints are found to be: 

1^24 + 1^34 _ 1^23 _ ^32 
UJ42 + 1^43 1^43 '-'^42 

provided that the height of the shocks is: 



(8) 



P 



UJ42 + (^43 



UJ42 + "^43 + "^24 + 1^34 



One should note that these constraints are exactly those obtained in [8] for the model 
with periodic boundary conditions. Under these constraints the time evolution equations 
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for \Pm,n) take the following form: 

H\Pm,n) = ^lr\Pm+l,n) + ^ll\Pm-l,n) + ^2r\Pm,n+l) + ^2l\Pm,n-l) — 

[Sir + Sii + 62r + S2i)\Pm,n) for m = 1, ■ ■ ■ , L - 2 and 
n = m + 2,- ■ ■ ,L 

H\Po,n) = -S\Pl,n) + 52r|Po,n+l) + (^2/ |Po,n-l) " 

(-5 + S2r + S2i)\Po,n) for n = 2, ■ ■ ■ , L 

(9) 

-f^|-Pm,L+l) = Slr\Pm+l,L+l) + l^li |-Pm-l,L+l) + S\Pm,L) — 

{Sir + Sii + 6)\Pra,L+i) for m = 1, ■ ■ ■ , L - 1 

H\Pm,m+i)=0 for m = 0, ■■■,L 
in which the bulk hopping rates for the shock fronts (random walkers) are: 



P ' 

il-p? 

and also the reflection rate from the boundaries is found to be: 

5^ = 52.-^^52.. (11) 

Oil + 02r 

Note that the original model has six independent parameters while by imposing the 
constraints ([H]) only four independent parameters are remained. 

The last equation in ([9]) implies that an empty lattice is a trivial steady-state of the 
model. It should be mentioned that the random walkers repel each other so that their 
positions cannot get closer than a single lattice site. This has already been observed in 
[7j for the special choice of parameters given in ([6]). For the present case, as in |7j, one 
can simply define a new measure \Pm,n) as: 

o) nO no) (12) 

for < ?Ti < 72 — 1 and 1 < n < L + 1 

and investigate its dynamics to see that the shock positions never meet each other at 
two consecutive sites. The reason for this goes back to the fact that the dynamical rules 
([5]) do not allow an empty system to be generated from a system with particles. In fact 
the system has two different steady-states. The trivial steady-state, as we mentioned 
above, is the one without any particles. The other one, which is a nontrivial steady- 
state, contains some particles. If we start with an empty lattice, it always remains 
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empty. In contrast, if we start with a partially- filled lattice, it will reach to an steady- 
state which contains particles. It can be shown that the nontrivial steady-state of the 
system can be written as a linear superposition of |Pm,n)'s. This will be discussed and 
presented elsewhere. In the next section we are going to calculate the steady-state of 
the system using the matrix-product approach. The question is that whether or not a 
finite-dimensional matrix representation exists for the quadratic algebra of this model 
under the constraints ([8]). 

4. Matrix-product steady-state 

According to the matrix product formalism the stationary probability distribution P{r) 
of any configuration t = {Ti\i = 1, ■■■, L} is assumed to be of the form: 

P{r) = ^{W\ f[{r,D + (1 - n)E)\V) (13) 

in which the occupation number Tj is defined as follows: Tj = if the site i is empty 
and Tj = 1 if it is occupied by a particle. The factor Zl in flT^ is a normalization factor 
and called the partition function of the system and can be easily calculated using the 
normalization condition J2{t,} P{'^) = 1- For a system described by the Hamiltonian 
of type ([2]), the operators D and E, which stand for the presence of a particle and an 
empty site at each lattice site respectively, beside the vectors {W\ and \V) should satisfy 
the following quadratic algebra: 




IWMl 



<W\\ - \=Q, [ \\V >=Q 






(14) 



in which the operators E and D are two auxiliary operators and h is of the form ([3]). 
Using flT^ one can simply calculate the corresponding quadratic algebra of our model: 

[C,C] = [E,E]=0, 

(Cc'32 + CJ42 + UJ2i)EC - (CJ32 + UJ42 + ^0-24 - ^23)-£^^ " 0J24C'^ + 

(cu24 - uj2z)CE = EC~ EC, 

(^34 - ^32)^C' - (u;23 + ^43 + ^34 " ^^32)^'^ - ^^34^^+ (15) 

(CJ34 + CJ23 + ^iz)CE = CE- CE, 
{W\E= {W\C = 0, 

E\V) = C\V) = 

in which we have defined C := D + E and C := D + E. 

For the special tuning of the parameters given in ([6]) the authors in [19] have shown. 
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for the first time, that the quadratic algebra (TT3]) has a four-dimensional matrix 
representation with non-diagonal E and C. In our generalized model it turns out that 
the quadratic algebra ( TT5l) also has a four-dimensional matrix representation provided 
that the constraints ([8]) are fulfilled. Note that these constraints are automatically 
fulfilled for (E]). 

As long as the eigenvalues of the matrix C are not equal and in an appropriate basis 
this matrix can be brought into a diagonal form. We have found that the matrices 
C and E besides the auxiliary matrices C and E have the following four- dimensional 
representation: 



C 



c 



E 



/ Sir <52r 

<5ii S21 




V 
1 ° 


(1- 




-pyt 







\ 


S2r 

021 

l) 



5 

o\ 




p^ 


-('^1. 


-(1- 


p)5ii) 















; 


VO 

<5i; &21 




(1- 




ei2 
-Pft 


0/ 

ei3 ei4 ^ 

624 










(1-P)|7 ^34 


V 
/ 






ei; 


1 ) 

5 ei4 \ 
















-P|7(^i^ - 


{l-p)5ii) 634 


V 









( 


3 ) 



(16) 



E 



in which 612, Cia, 613, eu and 634 are functions of cujj's and are presented in the Appendix 
A. In this representation we have chosen the other parameters 614, 624 and 634 in (TTBII as 
free parameters. On the other hand, the vectors \V) and (yV\ have the following matrix 
representations: 



1 

V^4 / 

in which Vi and 1^4 are free parameters, while v^ and 1^3 are complicated functions of 



1^) 



(iy| = (1,0,^3,^4 



:i7) 



UJ. 



1-J 



's and are presented in the Appendix A. Nevertheless, if one looks for a nontrivial 
steady-state then the relation between vi and W4 should be obtained by requiring that 
the probability of finding an empty lattice in long time limit is zero i.e.: 

{W\E^\V)=0. (18) 
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Figure 2. The phase diagram of the model in terms of the two variables x and y. The 
four regions I-IV are defined in the text. The coexistence line between the regions II 
and III is y = i. 



By considering the constraint flTSj) the partition function of the model can be easily 
calculated as follow: 



, 5l T. 5' 



Ir <^'2r^L 



2^2r.L 



a,C-i^:^r + a^iil - pY-^ 



,5. 



021 



2r^L 



+ 04 



(19) 



1^23 


^ ^32 


UJ2A 


^ C^34 


UJi2 


^^ ^43 


i — 


^L-i 



hi hi hi 

in which the coefficients a^'s are function of tUjj's. The explicit form of these coefficients 
are presented in the Appendix B. The phase diagram of the model can now be obtained 
by analyzing the thermodynamical behavior of the partition function Zl (i.e. in the 
limit L -^ oo) which obviously depends on the values of the hopping rates of the shock 
positions i.e. 5ii, (5ir, 821 and 52r- One should note that the model has a mirror symmetry 
which means that it is invariant under the following transformations: 



(20) 



Defining two new variables x := j^ and y := -f^, it can be shown using (TTOj) that the 
case {x > l,y < 1) never happens because p then has to be negative. In terms of the 
two variables x and y the phase diagram of the model has four regions defined as: 

Region I for x > l,y > l,xy > 1, 
Region II for x < l,y > 1, xy > 1, 
Region III for x < l,y > 1, xy < 1, 
Region IV for x < l,y < 1, xy < 1. 

The phase diagram of the model is shown in Figure 2. One can easily check that under 
the symmetry fl20l) the two variables x and y become 1/y and 1/x respectively. This 
implies that we need only to study two phases (either I and II or III and IV) since the 
information for the two other phases can be easily obtained using (!20l) . In what follows 
we only consider the two phases III and IV. In the phase III the left random walker (the 
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left shock front) has more tendency to move leftwards while the right random walker 
(the right shock front) has more tendency to move rightwards. In this phase the mean 
distance between the two random walkers is of order L and therefore the density of 
particles in the bulk of the lattice increases (this can be understood by looking at the 
definition of \Pm,n))- We call this phase the high-density phase. In contrast, in the 
phase IV both random walkers have more tendency to move leftwards, and therefore, 
the mean distance between the random walkers (the shock fronts) decreases. In this 
phase the mean distance between the two random walkers is of order 1. Since the mean 
density of particles decreases in this phase, it is called the low-density phase. Since the 
matrix representation of the quadratic algebra is finite-dimensional, the density profile 
of the particles defined as: 



ipk) 



for l<k<L (21) 



Zl Zl 

either has linear or exponential behaviors. Note that in calculating (12T1) the constraint 
flTSj) should also be applied to its nominator. It turns out that fl2T|) has only exponential 
behaviors in both phases III and IV. In the low-density (high-density) phase the mean 
density of the particles in the bulk of the lattice is zero (p) while at the boundaries 
the density profile of the particles changes exponentially with two correlation lengths 
I lny|~^ and | lnxy\~^ (| lny|~-'^ and | Inxl^""^). 

It is also interesting to study the coexistence line between the two phases III and IV i.e. 

One should note that in this case two eigenvalues of the 



on the line 



hr 



Hot ^ < 1. 



matrix C become equal, and thereby cannot be diagonalized. On the coexistence line 
between the two phases III and IV, we have found the following matrix representation 
for the operators in (ITBl) : 



C 



C 



E 



E 



( t 


\ 






(1-p) 










1 C34 


•) 


V 


1 ) 




/ 


^ 


\ 


Pi6,r-il 


-p)6ii) 










•) 


[o 


0, 


1 


(t -- 


ei3 ei4 \ 




(l-p)2 


624 







(1-p) 634 


5 


V 


1 ) 




/ 


ei3 ei4 \ 








-p{5,r 


-{l-p)5ii) 634 


V 







) 



(22) 
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and for the vectors: 



\V) 



{W\ = {1,0, ws,w^) (23) 





1 

which is vahd if the reaction rates uji/s in ([5]) satisfy both ([8]) and the following 
constraint: 

1^43 = UJ42[ ) + l^34( ) • (24) 

1-1^42 + '^43 + ^^24 + '^34 "^42 + "^43 + "^24 + "^34 

As in the diagonal case some of the matrix elements are free. In fl22l) the free parameters 
are chosen to be C34, 614, 624 and 634. The rest of the matrix elements are functions of cuj/s 
and are presented in the Appendix A. The two elements vi and w^ can be eliminated 
from the calculations by applying the constraints (fTSll while the other elements V4 and 
^3 are presented in the Appendix A. The partition function of the system, after applying 
the constraint flTSl) . is obtained to be: 

Zl = h{^f + 62(1 - pf- + hL + 64 (25) 

dii 

where the coefficients 6j's are complicated functions of the reaction rates cUjj's. These 

coefficients are presented in the Appendix B. As can be seen, on the coexistence line 

the left random walker (the left shock front) has more tendency to move leftwards while 

the right random walker (the right shock front) performs an unbiased random walk. In 

this case we expect that the density of particles increases near the left boundary. The 

density profile of the particles, defined in (12T|) . can be calculated on this line and it 

turns out that it is an exponentially increasing function near the left boundary with a 

characteristic length | lna;|~^ while it changes linearly in the bulk of the lattice. 

5. Summary and outlook 

In this paper we considered a generalized coagulation-decoagulation model on a finite 
lattice with reflecting boundaries and studied the conditions under which the steady- 
state probability distribution of the system can be written as a product of non- 
commutative operators. We started with a product shock measure with two shock fronts 
and obtained its dynamics generated by the Hamiltonian of the system. We found the 
necessary conditions on the microscopic reaction rates under which each of the shock 
fronts has a simple biased random walk dynamics. Finally by studying the steady-state 
of this system using the matrix-product approach we showed that its quadratic algebra 
has a four- dimensional matrix representation under the same conditions. The phase 
diagram of the model and also the density profile of the particles in each phase and on 
the coexistence line were studied. We saw that the whole phase diagram of the model 
can be described by the hopping rates of the shock fronts. The system has three different 
length scales. These are given by the eigenvalues of the matrix C which can be written 
in terms of the hopping rates of the shock fronts. 



Finite- dimensional representation of the quadratic algebra 12 

One can easily check that for the special tuning of the parameters (jHD the constraints 
(IH]) are satisfied and that all of our calculations reduce to those obtained in [19]. It 
is interesting to investigate the solutions of ([9]) in long time limit and construct the 
steady-state of the system as a hnear superposition of such solutions and then compare 
the results with those obtained from the matrix-product approach in present work. The 
results will be presented elsewhere. 

Appendix A. The matrix elements 

Here we present the matrix elements for the diagonal matrices (TT6l) and the vectors (TT7|1 . 
The matrix elements of E are obtained to be: 

ei2 = ((ei4(^24 + ^^34) ^^43 (^24 (^42 " ^43) + tu|4(^42 " ^43) + 2^^34(^42 (CJ42 + ^^43) + 

^^42(^42 + ^43)^ + ^24(C^42 - UJ43){2uJu + UJ42 + ^^43)) (^14(^42 " ^^43) + ^|4(^42 - ^^43)" 

^^43 (^42 +^^^43)^ + ^34(^^42 - ^43) + 2u;24(t^34(^42 -^^43)" 

'^43('^42 + ^43))))/(e24(t^|4^42 + 2u;24t^34t^42 + ^34.(^34(^^2 + 

((JJ42 + UJ43y)){^42{'^24{-UJ42 + ^^43) + UJ43{2uJ34 + UJ42 + ^^43 ) ) (CU24CU42 + 

2a;24t^34t^42 + tU34(tf 341:^42 + (^42 + ^^43)^)) + ^43(^42(2a;24 + ^34 + W42) + ("^34 + 

'^42)a;43)(t^l4^43 + ^^14^^43 + t^24 ( 2a;34a;43 + (^42 + Ci;43)^))))), 

ei3 = ((ei4CU43(u;24 + C^34 + UJ42 + C^43)^ (^14(^42 " ^^43) + ^34(^^34(^42 " ^43) + 
tc;42(w42 + ^43)) + C^24(2u;34(c^42 " ^43) " ^^43(^42 + ^43)))^)/(e34(c^42 + ^^43) (cu|4U;42 + 
2a;24CU34t^42 + tU34(t^34t^42 + (^42 + ^43)^)) (^14^42(^42 " 2^43) + t^34t^43(^|4(-2c^42 + 
a;43) - 3u;34U;42(c^42 + ^^43) - ^^42(^42 + ^^43)^) + U;i4(-3u;42^43(t^42 + UJ43) + 
U;S4{2ujl2 - 6U!42U}43 + uj^)) + t^24( "^42^43(^^42 + ^^43)^ + 

^34(^42 + ^43)(i^l2 - 4u;42u;43 + ujI^) + u;|4(cu|2 - 6a;42^43 + 2ujI^))))). 
The other matrix elements ei4, 624 and 634 are free. For the matrix E we have found: 
ei3 = ei3u;42, 

ei4 = -(ei4U;42Cc'43(^24 + t^34 + t^42 + t^43)^(^42((^24 + ^^34)^ + ^34^42)- 

((u;24 + ^^34)^ + (u;24 - ^^34)^42)^^43 - C^24t^|3))/((^42 + ^^43) (t^24^|2((^24 + ^34)^ + 

^^34^^42) - u;42(2(a;24 + ^^34)^ + 3(u;|4 + u;24^34 + ^14)^42 + (^^24 + ^^34)^^2)^43+ 

(^34(^^24 + C^34)^ - 3(u;|4 + U;24t^34 + ^^14) ^42 " 2(u;24 + ^^34)^12) ^43 + (^24^^34" 
(Cc'24 + ^34)^^42)^43))) 



634 — — 634^43. 
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For the vector elements in (1171) we have found: 

V4 = (t^i4(c^42 - UJ43) + ^^34(^^34(^^42 - ^^43) + C^42(C^42 + ^^43)) + 
^24(2u;34(^42 - ^43) - ^43(^^42 + ^43)))/(e34(w24^42 + 26^24^34^42 + 
U;34(u;34Ct'42 + (C^42 + a;43)^))), 

^^3 = -(ei4^42(t^24 + ^34 + ^42 + ^43)^(t^42((c^24 + ^^34)^ + U;34U;42)- 

((t^24 + ^34)^ + (C^24 - ^34)^^42)^^43 - ^24^13)) 7(^34 (^42 + ^^43) (t^24C^|2((^24 + ^34)^ + 

^34^42) - Co'42(2(t^24 + ^0-34)^ + 3(u;|4 + UJ24UJ34 + ^34)^^42 + (^24 + ^34)^|2)^43 + 

(^^34(1^24 + ^^34)^ - 3(60-^4 + 60-24^34 + 0ol4)UJ42 - 2(60-24 + ^34)^^12)^43 + 

(t^24t^34 - (^24 + ^34)^^42)^43))- 

On the other hand for the non-diagonal matrices (122!) and the vectors (123!) we have 
found the following results. The matrix elements of E are obtained to be: 

ei2 = (-£34614(^12 - ^^13) (^34(^42 - ^43) ^ + (^42 - ^43) (-^34^12 - 6a;34a;42C^43- 

"^43^34) + ^43^42 + "^421^43 + ^^42^43, + Yi^43(l-^34 + 1^42) " 1^341^42 ( — 2ti;34(u;42 — ^^43) + 
4^^43^42(^42 + ^43))))/(e24(\/^43(w34 + a;42) " ^^34^^42 (634^43 (^|2 (^34 " 10u;34^42 + 
5^42) - 2(t^34 - 5u;42)cc'42(t^34 + ^^42)^^43 + (^34 + ^42)^^13) + C34W42 (^34 (^42 -^^43)^" 
2a;34(cJ42 - CU43)(3u;42 + U;42^43 + 2u;43) + ^42(^12 + 6u;|2^43 + 5UJ42UJI3 + 4:UjI^))) + 

e34u;42u;43(t^l2(5ct;i4 - 10u;34C^42 + UJI2) + 10u;42(-u;|4 + ujl2)^i3 + 5(^34 + ^42)^^13) + 

£34(^34(^42 - ^43)^(4^^42 + ^43) - 2u;34U;42(u;42 - UJ43){^42 + CJ43)(2u;^2 + 2u;42u;43 + 



UJ 



43, 



+ ul2UJ43{4:Ujl2 + 5co'42Co'43 + 6uJ42UjI^ + ^Is)))), 



ei3 = (-ei4a;43(w42 - ^43))/((c34^42 + 634^43)^/^43(^^34 + C^42) " ^34^42 + 
e34'^43'^42 + C34Ci;43). 

As in the diagonal case the other matrix elements ei4, 624 and 634 are free. The matrix 
elements for E are found to be: 

^13 = ei3a;42, 

ei4 = ( (6146346^42 (u;34 + Cc'42)t^43))/(C34(^43 + ^^42) y^t<;43(u;34 + UJ 42) - i034L042^ 
634(^34 + ^42)^43 + 634^^42(^^34 - ^43)), 

634 = — 634a'43. 

The non-free vector elements in this case are given by: 

V4 = {uj42 - ^43)7(634(^^42 + y ^43(^^34 + ^42) - ^34^42)) 

W3 = (6i4u;42(u;34 + u;42)) 7(634^43(^34 + ^^42) + 634(^42(^^34 - ^43) + 



(cc;42 + u;43)^u;43(u;34 + u;42) - ^34^42)) 
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Appendix B. The partition function 

The coefficients a^'s for z = 1, ■ ■ ■ , 4 in (iT9l) are found to be: 

ai = -{ei4UJ42{uJ42{{uJ24 + ^34^ + ^^34^^42) - ((^^24 + ^^34)^ + (^24 " ^34)^^42)^43- 
^24^43) (^42 ((t^24 + ^^34)^ + ^^34^^42) - ((t^24 + ^34)^ + 2u;24C^42 + ^^42)^43 " (2c^24+ 
t^34 + 2UJ42)UJJ^ - ^43) (^24 + ^34(^34 + ^42) + ^24(2(u;34 + ^42) + tU43)))/(e34(^42 + 
^43) (^42 (-(^24 + ^34)^ + ^24^42) + ((C^24 + ^^34)^ + 2a;24C^42)^43 + ^24^43)1^24^^42 
((a;24 + "^34)^ + '^34'^42) - Co'42(2(co'24 + ^^34)^ + 3(u;|4 + 0^24^34 + ^|4)^42 + (^24 + 
a;34)a;|2)c^43 + (^34(^^24 + ^^34)^ - ^{^24 + ^24^^34 + tu|4)^42 - 2(cJ24 + ^^34)^42)^43 
+ (t^24t^34 - (^24 + ^^34)1^42)^13)), 

^2 = -(ei4(c^24 + t^34)^t^42('^42(('^24 + ^^34)^ + ^34^42) - ((^^24 + ^^34)^ + 

(u;24 - ^^34)^42)^43 - UJ24UjJs){uJ42{uJ24 + (^34 + ^42)^ + ^24(2(^34 + ^^42))" 

((u;24 + ^^34)^ - 2u;34U;42 - 2u;^2)^43 - (^24 - ^^42)^^13) (^42 ((t^24 + ^^34)^ + ^^34^^42)" 

((u;24 + C^34)^ + 2t^24t^42 + ^^42)^43 " (2t^24 + ^^34 + 2cJ42)u;|3 - Cc'|3))/(e34(c^42 + a;43) 

(^42(-(^24 + ^^34)^ + ^^24^^42) + ((^24 + ^34)^ + 2^^24^42)^43 + ^24^43) 

(^42 ((^24 + ^^34)^ + UJ34UJ42) - {{UJ24 + ^1^34)^ - 2u;34U;42)^43 + ^34^^13) 

(u;24^|2(-(^24 + ^34)^ - ^34^^42) + ^42(2(u;24 + ^34)^ + 

3(u;|4 + U;24^34 + ^34)^42 + (^24 + ^^34)^12)^43 + (-^34(^^24 + C^34)^ + 

3(co'|4 + U;24^34 + ^34)to'42 + 2(cc'24 + ^34)^^42)^43 + ("^24^34 + (^^24 + ^^34 ) ^42 ) ^^43 ) ) , 

^3 = -(ei4^42(^24 + ^34 + ^42 + ^43)^(^42((t^24 + ^^34)^ + ^^34^^42) " ((^24 + ^34)^ + 
{UJ24 - CU34)u;42)^43 " tU24t^|3))/(e34(^42 + ^43) (^24^12 ((^24 + ^^34)^ + ^^34^^42)" 
'^42(2(CU24 + C^34)^ + ^{^h + ^24^34 + ^14)^42 + (^24 + ^^34)^12)^43 + (^34(^24 + 
^34)^ - 3(^0-14 + U;24t^34 + to'34)'^42 " 2(u;24 + "^34)^12)^43 + (^24^^34 " (^24 + 



04 = (ei4U;42(cc'42((^24 + ^^34)^ + ^34^^42) " ((^^24 + ^^34)^ + (^24 " ^^34 ) ^42 ) ^43 " 

1^241^13) i^'i2{^24 + (^34 + ^42)^ + ^24(2u;34 + UJ42)) - ((^24 + ^34)^ - 2u;34U;42- 

2a;|2)'^43 - (^24 - ^^42)^^13) (^24 + ^24(2^^34 + ^43) + ^34(^^34 + tU42 + 

2a;43)))/ (634(^42 + W43)('^42((tt'24 + ^^34)^ + ^34^42) " ((^^24 + ^^34)^" 

2a;34a;42)c^43 + ^34^^43) ('^24a;|2( ('^24 + ^34)^ + ^^34^^42) - ^42(2(u;24 + ^34)^ + 

3(u;|4 + CJ24^34 + 1^14)^42 + (^24 + ^^34)^42)^43 + (^^34(^24 + ^34)^- 

3(u;|4 + U;24^34 + ^14)^42 - 2(u;24 + ^34)^^12)^43 + (^24^34 - (^^24 + Ct;34)a;42)t^|3)), 

The coefficients 6,'s for i = 1, ■ ■ ■ , 4 in (l25l) are also found to be: 
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h = (-034614^^42^43(^12^43(^42 + ^43)(3CJ42 + ^43)^ + ^14(^42 - t^43)^(8c^|2 + 
^^42^^43 + OOJ^) - 2CU34CU42(CJ42 " ^43)(4a;|2 + II0JI2UJ4S + 4^42001^ + ^l3) + 

Cc;43(c^34 + 6^42) - ^^34^42(26^42(^^34(^^42 - 6^43)^ + 2CJ34(-3CJ42 + Cj|2^43 + ^^42^^43 + 



^h) + C^42(tc|2 + 8a;|2^43 + 5CJ42CJ|3 + 2ujI^))))) / {€34(^12 - C^|3) ((-^34^^42 + (^34 + 
^42)^43)(4e34U;42C^43(^34(-C^42 + 6^4 3) + ^42(^^42 +^^43))+ £3 4(^42(^42 + ^43)^ 
-C^34(C^42 - C^43)(3u;42 + Ujjs))) + ^^isi^Si + ^^^42) " ^^34^^42 (634^^43 (^^42(^34" 

6a;34a;42 + UJI2) - 2(cj34 - 3CJ42) 6^42(^^34 + ^^42)^^43 + (^34 + (^42)'^ (^43) + 

C34UJ42{l^u{^A2 - ^^43)^ - t^34(tc'42 " ^0^43) (36^42 + 2cJ42U;43 + 3c^43) + 
CJ42^43(3u;|2 + 2CJ42C^43 + 3CJ43))))), 

h = (£34614^^42(^42 + ^43) (^34(^42 - t^43)^(3c^42 + 2CJ4 3) + 6^126^43(^42 + 10C^42^ 43 + 
bujj^) - U;34C^42(C^42 " ^43) (^|2 + I2CJ42CJ43 + 7ujI^) + ^^(^43(^34 + ^^42) " ^34^^42 (^^14 
(c^42 - 6^43)^ - C^34(c^42 - C^43)(3cj|2 + 8u;42CJ43 + CJ43) + iO 42UJ 43{hujl2 + 10u;42CJ43 + 
'^43))))/(e34(t^42 - ^^43) ((-^34^^42 + (^34 + ^^42 ) ^43 ) ( " C34 ( 5cJ34 " CJ42)u;|2 + 
^42(634(-7cJ34 + 5c^42) + C34(4c^34 + 6a'42))c^43 + C^42(-(c34 " 6634)^^34 + 

5 (c34 + 2634 )a;42) 0^43 + (63 4(^^34 + 6^42) + 2c34(a;34 + 2uj42))ujI^) + 

\J^iZ\!^M + 1^42) — '^34'^42 (634^43 (2^34(^^42 ~ ^43) — 
C^34U;42(C^42 - C^43)(9C^42 + 7u;43) + ^|2(^42 + 
10tc;42C^43 + 5a;43)) + C34(2ci;|4Ci;42(cc'42 - c^43)^- 
^34(^^42 - ^43)(4CJ42 + 6cj|2C^43 + 5cJ42C^43 + (^43) + 
^42^43 (4C^|2 + 5C^|2C^43 + QUJ42UJI3 + 6^43))))), 



^3 — \c^A^iA^42\!^42 — ^43))/ (634(634^43 + 6340^42 ) y 0^43 ((^34 + 0742) — u;34a;42+ 

e34U;43 (634^^42 + e34U;43)), 

64 = {—C34ei4UJ42{uj42 — (^43) (4^^42(^341^42 " (l^34 + 1^42)^^43) (^34 (t<-'42 — ^43) — 
t^42(t^42 +^^43)) + Y ^43(^^34 + C^42) - ^34^42 (a;|2('^34 " 6CJ34CJ42 + CJI2)- 
2(a;34 - 3cU42)tc'42(u;34 + ^^42)^^43 + (^34 + t^42)^Cu|3)))/(e34(a;42 + U;43)((3CJ42 + CU43) 
(-U;34a;42 + (0^34 + C^42)t^43)(e34t^42(t^42(-t^34 + ^^42) + (t^34 + ^^42)^^43 + 2c^|3) + 
634^^43(^34 (-^42 +^^43)+ C^42(3u;42 + U;43))) + (- 26^34(0^42 -1^43) + 
U}42{UJ42 + 3c^43))(634U;42C^43(t^42(-3u;34 + ^^42) + 3(cJ34 + U;42)c^43) + 
e3 4(-U;34(c^42 - C^43)(2u;^2 + ^43) + ^42^^43 (2C^42 + ^^42^^43 + C^|3))) + 
Cc^43(t^34 + ^^42) - ^34^^42 (26346^43 (0^34(0^42 - 6^43)^ + ^^42(3^42 + 



C^43)(^42 + 3CJ43) + 4a;34C^42(-2c^|2 + 6^42^43 + ^43)) + 

C34(2a;|4C^42(c^42 - ^43)^ + ^34(-9c^42 + 2u;|2u;43+ 

QUJ42^43 + ^43) + ^42(^^42 + 10u;|2CJ43 + 
100^12^^43 + 10^42^43 + ^43)))))- 
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